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ABSTRACT 



A multiple carrier communication system includes a primary 
impulse shortening filter that receives an output signal of an 
analog to digital converter and accepts coefficients. A sec- 
ondary impulse shortening filter receives the output signal of 
the analog to digital converter, outputs an output signal, and 
passes coefficients to the primary impulse shortening filter. 
A reference signal generator outputs a reference signal. A 
comparator compares the: output signal and the reference 
signal and outputs a resulting error signal. An adaptive 
processor computes coefficients for the secondary impulse 
shortening filter based on the error signal. 

28 Claims, 12 Drawing Sheets 
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TIME DOMAIN EQUALIZATION FOR 
DISCRETE MULTI-TONE SYSTEMS 

CROSS REFERENCE TO RELATED 
APPLICATIONS 

This application claims priority from U.S. Provisional 
Application No. 60/087,336, filedMay 29, 1998 and entitled 
"Time Domain Equalization in DMT Systems", which is 
incorporated herein by reference. This application claims 
priority from U.S. application Ser. No. 09/054,468, filed 
Apr. 3, 1998 and entitled "Spectrally Constrained Impulse 
Shortening Filter for a Discrete Multi-tone Receiver", which 
is incorporated herein by reference. 

BACKGROUND 

The invention relates to time -domain equalization in a 
discrete multi-tone (DMT) receiver. 

Conventional single carrier modulation techniques trans- 
late data bits for transmission through a communication 
channel by varying the amplitude and/or phase of a single 
sinusoidal carrier. By contrast, DMT, which is also referred 
to as Orthogonal Frequency Division Multiplexing (OFDM) 
or Multicarrier Modulation (MCM), employs a large number 
of sinusoidal subcarriers, e.g., 128 or 256 subcarriers. The 
available bandwidth of the communication channel is 
divided into subchannels and each subchannel communi- 
cates a part of the data. A DMT system may employ 
quadrature amplitude modulation (QAM) for each of the 
subcarriers. 

OFDM-based systems transmit blocks of information bits. 
The time required to transmit one such block is called the 
symbol period. The time domain waveform that corresponds 
to one such block of bits is called a symbol. 

Intersymbol interference (ISI) arises from the character- 
istics of practical communication channels and limits the 
rate at which information can be transmitted through them. 
Specifically, communication channels typically have an 
Effective Discrete-Time Impulse Response (EDIR) that is 
greater than one sample time in length, which causes ISI. ISI 
is a well-known phenomenon in single-carrier communica- 
tion systems and there are many techniques for reducing it. 
The process of such ISI reduction is called equalization. ISI 
is discussed, for example, in Proakis, Digital 
Communications, McGraw Hill, 2nd Edition, 1989. 

Equalization in OFDM-based systems is achieved by a 
two stage process. First, at the transmitter, a Cyclic Prefix 
(CP) is employed by affixing an end-portion of each symbol 
to the beginning of the symbol. A cyclic prefix that is greater 
than the EDIR of the channel prevents one symbol from 
interfering with another. Furthermore, it also facilitates a 
simple method of neutralizing the time domain spread of 
each symbol forced by the channel. This is achieved through 
a simple frequency domain process in the receiver which 
requires one multiplication operation for each subcarrier 
used. The use of a Cyclic Prefix to reduce ISI is discussed, 
for example, in: Cimini, "Analysis and Simulation of a 
Digital Mobile Channel using Orthogonal Frequency Divi- 
sion Multiplexing," IEEE Transactions on communications, 
pp 665-675, July 1985; Chow, "A Discrete Multi-Tone 
Transceiver System for HDSL applications," IEEE Journal 
on Selected Areas of Communications, 9(6): 895-908, 
August 1991; and "DMT Group VDSLPMD Draft Standard 
Proposal/' Technical Report, T1E1.4/96-329R2, ANSI 
1997. 

Another problem arising in conventional DMT systems is 
noise bleeding, which occurs when noise from one fre- 
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quency band interferes with a signal whose subcarrier is in 
another frequency band. Noise bleeding is caused, in 
general, by a discrete Fourier transform (DFT) operation at 
the receiver. Noise bleeding is discussed in, for example, 

S Worthen et. al, "Simulation of VDSL Test Loops," Techni- 
cal Report T1E1 .4/97-288, ANSI 1997. 

In a perfectly synchronized DMT system, a signal in one 
frequency band does not interfere with a signal whose 
subcarrier Is in another frequency band. However, noise 

10 from one band may interfere with other less noisy bands and 
render them unusable. Techniques for dealing with noise- 
bleeding include wavelet-based solutions. However, 
wavelet-based solutions are, in general, computationally 
intensive. 

Other references dealing with time domain equalization 
include: Chow, J. S. and Qoffi, J. M., "A Cost-effective 
Maximum Likelihood Receiver for Multicarrier Systems", 
Proceedings of the ICC, 1992; Melsa, Peter J. W., Younce, 
Richard C, and Rohrs, Charles E., "Optimal Impulse 
Response Shortening"', Proceedings of the thirty -third 
20 Annual Allerton Conference on Communication, Control 
and Computing, 1995, pp. 431-438; Harikumar, Gopal and 
Marchok, Daniel, "Shortening the Channel Impulse 
Response of VDSL Loops for Multicarrier Applications", 
Technical report T1E1 .4/97-289, ANSI, 1997. 

25 SUMMARY 

^^f^ll^^^^j^pl^ffi^in^^^SM^y 

jipjjtKe lSjll^ 
30 c|geffi'cielTts^^ 
mg^io^ss:;thlffl 
c cjtion/systenKM 
(gfeSMSFi^ 
/jritcrs^f^j. Interfere 
35 eifective'dLscr^t'e^tirhe impulse response^HDIRYoTthe com- 
^ munication channel: Conventional-impulse shortening 
ma^ have"deepTnulb~m 

contrast, the' SClSF"has"a filter characteristicr that is" essenr 
^allyjff^^ 

40 p^teljr^Hmmite^ 

Second, the SCISF reduces noise bleeding between sub- 
channels by attenuating noisy channels in a manner that does 
not reduce the signal to noise ratio (SNR) in these channels, 
but reduces the noise power that may appear in the sidelobes 

45 of adjacent subchannels./^ ^^ gSfeaccomplishesrthese, 
fin^tiola^ 

p^edsbn^afcdesirej^pec tral response . 

The coefficients of the SCISF are computed independent 
of the length of the cyclic prefix, the symbol length, and the 

50 frequency domain equalization characteristics of the system. 
The SCISF is particularly effective in systems in which 
additive noise dominates intersymbol interference, or in 
which noise bleeding predominates. The SCISF reduces 
noise bleeding with a filter structure that is shorter than that 

55 obtained with other techniques. Consequently, the SCISF is 
less complex and its coefficients are easier to compute, 
which reduces system complexity and cost. The SCISF may 
be particularly well suited, for example, for very high-speed 
digital subscriber lines (VDSL) systems, which generally 

60 have low intersymbol interference and tend to suffer from 
noise bleeding. 

In addition, dynamic selection of a cyclic prefix (CP) that 
maximizes the data throughput for a communication channel 
having a particular noise profile is provided. Dynamic 

65 selection of the CP allows the communication system to 
adapt to changing noise conditions in the communication 
channel. 
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In one aspect, generally, a primary impulse shortening generator outputs a reference signal. A comparator compares 

filter is adapted in a multiple carrier communication system. the output signal and the reference signal and outputs a 

A secondary impulse shortening filter is provided. An output resulting error signal. An adaptive processor computes coef- 

signal of the secondary impulse shortening filter is compared ficients for the secondary impulse shortening filter based on 

to a reference signal to compute an error signal. Coefficients s the error signal. 

of the secondary impulse shortening filter are computed in Embodiments may include one or more of the following 

an adaptive processor based on the error signal. Coefficients features. A discrete Fourier transform may receive an output 

of the primary impulse shortening filter are replaced with s i gna i 0 f the primary impulse shortening filter. A decoder 

coefficients of the secondary impulse shortening filter. may reC eive the transformed output signal from the discrete 

Embodiments may include one or more of the following 10 Fourier transform. The reference signal generator may 

features. An output signal of the primary impulse shortening include an encoder that receives output data from the 

filter may be decoded to form output data. The output data decoder. The reference signal generator may include a 

may be encoded to form the reference signal. A discrete scaling filter that scales the output data from the encoder 

Fourier transform may be applied to the output signal of the using a set of scaling factors. An inverse discrete Fourier 

primary impulse shortening filter prior to decoding the 15 transform may receive the scaled output signal from the 

output signal. An inverse discrete Fourier transform may be scaling filer, 

applied to the encoded output data in forming the reference [n aspect> a multiple carrier communication sys- 

signal. tem may include a spectrally constrained impulse shortening 

A digital signal may be received from an output of an filter that receives an output signal of an analog to digital 

analog to digital converter. The digital signal may be input 20 converter and accepts coefficients. A reference signal gen- 

to the primary impulse shortening filter. The digital signal erator outputs a reference signal. A comparator compares the 

may be delayed. The delayed digital signal may be input to output signal and the reference signal and outputs a resulting 

the secondary impulse shortening filter and the adaptive error signal. An adaptive processor computes coefficients for 

processor. the spectrally constrained impulse shortening filter based on 

The encoded output data may be scaled with a set of 25 the error signal, 

scalingjactonsjn^^ Embodiments may include one or more of the following 

£f acto rs may be _ determined by; m^ur^ features. A memory may store the reference . signal as a 

(power spectral den^y ^^ojn pu^ g_a^esir e d^spe^ tral^ predetermined signal. A discrete Fourier transform may 

^resrxjnsebasedonmemeasur^n^ise power, and conipu^ing> ^ receive the reference signal from the, memory. A scaling 
(the scaling factors so that the co^ffiaents computed in_the^ filter may scale the reference signal using a set of scaling 

(Jia^pu^pxpce^sor^prgvide the secondary ^uls^horte^ • factors. An inverse discrete Fourier transform may receive 
^g filteTwith a spectral response that mat^h^th"e~des^ifed - the scaled reference signal. 

spectral response. ^ATiicrete Fourier, transform^ may^be 0ther features ^ advantag e S will be apparent from the 

appi^ ^c 7 o«^nal^f the pri^y imr^ short, 3s followmg detailed description, mcMng the drawings, :and 

<ernng_fyter-prior-to-decoa^Jh^^ noise &Qm tte c \^ ms 
power spectral density may be measured at an output of the 

discrete Fourier transform. An inverse discrete Fourier trans- BRIEF DESCRIPTION OF THE DRAWINGS 
form may be applied to the scaled, encoded output data. 

InTffio.tfrerafj^^ FIG - 1 is a block dia e ra m of a.discrete multi-tone com- 

^m-a^lhple-camercom munication system with a spectrally constrained impulse 

txafiyxonstr ained. impulse. shortening.filter. An output signal shortening filter. 

of the spectrally constrained impulse shortening filter is FIG. 2 is a block diagram of a system for adapting a 

compared to a reference signal to compute an error signal. spectrally constrained impulse shortening filter in an initial 

Coefficients of the spectrally constrained impulse shortening 45 training process. 

filter are computed in an adaptive processor based on the FIGS. 3A and 3B are block diagrams of a system for 

error signal. adapting a spectrally constrained impulse shortening filter in 

Embodi ments ma yj nclud e one or more of the following an periodic or continuous training process, 

features. Thejeferencg_si)maj5 ^ FIG. 4 is a block diagram of a generalized adaptation 

stored^Xmemo ry in thecommunic ation s ystems A discrete 50 system that employs frequency scaling in the feedback loop. 

Fourier transformmay be applied to predetermined refer- FIG. 5 is a block diagram of a system for adapting a 

ence values to form transformed reference values. The spC ctraUy constrained impulse shortening filter in an initial 

transformed reference values may be scaled with a set of lraining process ±t% includes frequency scaling in ih& feed . 

scaling factors to form scaled. An inverse discrete Fourier j 00 p 

transform may be applied to the scaled values to form the « ^ T ^o , A , , n . , . r r 

reference sienal FIGS. 6A and 6B are block diagrams of a system for 

. , , r , adapting a spectrally constrained impulse shortening filter in 

Adata signal may bci received from lan output of an analog ^ riodic 0f contirmous trainiag proccss ^ includcs 

to digi al converter. The data signal may be input to the f scaJi in ^ feedback { 

spectrally constrained impulse shortening filter and the 2___ ' _ e 

„ JQ FIGS. 7A-7D are plots of simulation results for a discrete 

adaptive processor. 60 

T ... . , . . x . multi-tone system. 

In another aspect, a multiple earner communication sys- J 

tem includes a primary impulse shortening filter that FIGS, 8A-8D are plots of simulation results for a discrete 

receives an output signal of an analog to digital converter multi-tone system, 

and accepts coefficients. A secondary impulse shortening DESCRIPTION 

filter receives the output signal of the analog to digital 65 

converter, outputs an output signal, and passes coefficients to As shown in FIG. 1, a discrete multi-tone (DMT) com- 

the primary impulse shortening filter. A reference signal munication system 10 has a transmitter 12 and a receiver 14. 
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The transmitter 12 accepts an input data bit stream which SCISF filters the output {y k } of the receiver A/D 80. The 
passes through a constellation encoder 20. The encoder 20 coefficients are selected using an algorithm that minimizes 
divides the serial input bit stream into blocks of data. These the squared error between a reference sequence {u^} gen- 
blocks of data are further subdivided into smaller blocks crated by the receiver and the output of the SCISF {H k }. The 
corresponding to subchannels. Each of these smaller blocks 5 SCISF may be a finite impulse response (FIR) filter^ or an 
are used to compute a complex value representing a con- infinite impulse response (IIR) filter. The SCISF may be 
stellation point. Each constellation point corresponds to a trained following activation of the communication system or 
subsymbol. The subsymbols are then output by the encoder. periodically during operation of the system to compensate 
Taken together, the subsymbols constitute a symbol. for variations m me channel noise profile . 

Hie subsymbols are supplied to an inverse discrete Fou- 1fl . . . f A . f . 

rier transform (IDFT) 30, which may be implemented, for 10 iU ™ e m *? be P e ^™d using a variation of one of 

example, in a digital signal processor. Hie IDFT 30 outputs ^Lf* cal "^T? l\f^™> ™& f least mean squares 

N time samples of a symbol. The time samples are processed (LMS), normalized LMS, or recursive least squares (RLS). 

by a parallel to serial converter 40 to form a single stream For example, the following algorithm is a version of the 

of time samples normalized LMS algorithm in which b 0 , . . , , b^ and a 1? . , , , 

Following the parallel to serial converter 40, a prefix 15 V are me FIR and IIR parts of a SCISF having an impulse 

adder 50 adds a cyclic prefix to the beginning of each res P onse > & ™ e ^^n* G(z) is: 
symbol to reduce intersymbol interference (ISI). 

Alternatively, the cyclic prefix may be added in the parallel G[z) = bo^i^-^^ W 

to serial converter. After the cyclic prefix is added, the 2Q l-^r 1 a N zr* 
resulting signal passes through a digital-to-analog (D/A) 

converter 60 for transmission to the receiver 14 through a ™ j , r ^ • . jl-j^j-.l 

communication channel 70. An analog transmit filter 65 may , ,™ e .adaptation of coefficients a, and b.^is defined in the 

be included following the DIA converter to band limit the f ° llowm S f which a^) and b,<k) are the values 

transmitted signal coefficients during the Ir iteration. The parameter 

. ,i . ' . iL . * it _ * , 25 u is a predetermined constant with a value of 0.4, 

At the receiver 14, the signal passes through an analog- ^ r 

to-digital (AID) converter 80 and then through a spectrally 

constrained impulse shortening filter (SCISF) 90. A prefix ^ _ y b .( k ) yk . + y a mh ■> ^ 

stripper 100 strips the cyclic prefixes from the resulting ' tif ' *" ' j^i 

symbols and a serial to parallel converter 110 divides the 30 
stream of time samples into parallel signal paths that form 
the inputs to a discrete Fourier transform (DFT) 120. The 
DFT 120 converts the time samples into subsymbols. A 
frequency domain equalization filter 130 equalizes the sub- 
symbols. A decoder 140 converts the subsymbols into a data 35 
bits and outputs the resulting data. An analog receive filter 
75 may be included prior to the A/D converter in order to 
band limit the received signal 

As discussed above, a cyclic prefix is added to each 
symbol prior to transmission through the communication 40 
channel to reduce the effects of ISI. The cyclic prefix is 
added by copying the last v time samples from the end of a 
symbol and placing them at the beginning of the symbol. To 
eliminate ISI, the length of the cyclic prefix, v, is chosen to 
be longer than the effective discrete -time impulse response 45 
(EDIR) of the channel. However, because the cyclic prefix 
constitutes redundant data, increasing the length of the [ n a fi^t embodiment, the coefficients of the SCISF are 
cyclic prefix reduces the efficiency of the communication determined during an initial training period following the 
system. For example, in a system having N time samples per activation of the communication system using the LMS 
symbol and a cyclic prefix of v time samples, the efficiency 50 algorithm described above. A predetermined sequence of 
of the system will be reduced by a factor of N/(N+v). b i te Xjt j s mput ^ tne transmitter 12. The sequence of bits 
Efficiency may be maximized either by minimizing v or results in a sequence of real numbers {x k } at the input of the 
maximizing N. However, increasing N increases the q/a 60. As shown in FIGS. 1 and 2, the transmitted signal 
complexity, latency and computational requirements of the ^ filtered and noise-corrupted by the transmission channel 
system and at some point becomes impractical. Accordingly, 55 70, resulting in a received sequence {y k } at the output of the 
it is desirable to minimize v. j^q go in the receiver 14. The SCISF 90 filters and 

A spectrally constrained impulse shortening filter having transforms the received sequence {y k } into an output 
an impulse response, g(n), may be employed in the receiver sequence {x k }. The output sequence {it k } is compared (in a 
to minimize the length of the cyclic prefix by decreasing the signal comparator 205) to the predetermined sequence x^, 
EDIR of the effective communication channel, which 60 which is stored in memory 210 in the receiver. The corn- 
includes the transmit and receive filters, the impulse short- parison results in an error signal e fc that is input to the LMS 
ening filter, and the physical transmission channel. The use algorithm processor 215 along with sequence {x k }). 
of an impulse shortening filter is referred to as time domain ^The-traming_process-determines-coefficients~for— the^> 
equalization. Decreasing the EDIR allows a shorter cyclic 'SCISF so that the output {x J matches Jhe predetermined^ 
prefix to be used without increasing ISI. 55 sequence {x^} as closely as possible in a least squares sense, f 

The SCISF may be configured, by determining the filter r Le., the mean square error between^th~e^output-and=th*e> 
coefficients during a training or adaptation period. The (^edetermined-sequence-is^mmimized.. During Jhe jrainin g: 
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/process, the coefficients of the SCISF converge to values The output sequence {k k } of the secondary SCISF 300 is 

Ithat enable- the^SCISF-to reduce ISI~and additive noise. The compared (in a signal comparator 205) to the reference 

< resulting SCISF matches P a (w) in the frequency, domainin sequence {x k }, which is output by the encoding/modulation 

LAleast.squares-sense, where: cna in (320, 330, 340 and 350). As noted above, the received 

— > 5 sequence {y k } passes through a delay 305 before being input 

" p^io) = ^ H *^ f ( 3 ) to the secondary SCISF 300. The delay 305 compensates for 

l^. — — -^(w)|W(w)P"+5^(to)° the processing delay in the demodulation/decoding chain 

and the encoding/modulation chain. The comparison results 

— _ in an enor signal e* that is input to the LMS algorithm 

In the equation i above, S^-i^-^^pecM-density^ rocessor 215 . ^ ^ occss d6t6rmines coefficients 

the,mput ofUhe-transmU^Q/A^, ( § n ((«) is the ^ power £ * { } 

spectral density of the additive noise at the output of the A/D c 3 , ^ , , i. t - i 

80r and H(co) is the frequency response of the effective the reference < x *> 35 closel y possible in a leas 

discrete-time impulse response (EDIR) of the transmission squares ser^, i.e., the mean square error between me output 

'channel 70;transmit filter 65, and receive filter 75 measured and l ne reference sequence ^ rmnimized. Periodically, the 

between the input of the transmitter D/A 60 and the output 15 coefficients of the secondary SCISF 300 are copied to the 

bf_the^receiver A/DM. P rimar y SCISF 90 

Upon completion of the initial training of the SCISF, the Alternatively, as shown in FIG. 3B, the periodic training 

SCISF coefficients are fixed and the frequency domain may be performed with a single SCISF 90. In this 

equalizer (FEQ 130) of the receiver is trained using standard configuration, a received sequence {y k } is output by the A/D 

techniques for DMT receivers. Following the training of the 20 80 in the receiver 14. The SCISF 90 filters and transforms 

FEQ, the SCISF can be periodically trained or adapted the received sequence {yj into an output sequence {x'J. 

during operation of the communication system. Since it is The received sequence {y k } is also input to a delay 305. 

not efficient to repeatedly transmit a predetermined bit After the received sequence {y*} passes through the SCISF 

sequence during operation of the communication system, the 90, a data switch 360 is changed from position A to position 

periodic training process uses transmitted communication 25 B, allowing the delayed received sequence to make a second 

data to generate the reference and output sequences, as pass through the SCISF 90. An output switch 370 also may 

described below. be opened, so that data is not output during the training 

During operation of the communication system, a process. In addition, the SCISF coefficients are controlled by 

sequence of communication data bits x k is input to the the LMS algorithm during the training process, 

transmitter 12 (see FIG. 1). The sequence of data bits results 30 A reference sequence is computed as in the configuration 

in a sequence of real numbers {x*} at the input of the D/A of FIG. 3A. The data bits are input to an encoder 320, an 

60. As shown in FIGS. 1 and 3 A, the transmitted signal is IDFT 330, a parallel to serial converter 340, and a prefix 

filtered and noise-corrupted by the transmission channel 70, adder 350. The resulting reference sequence {x k } is input to 

resulting in a received sequence {y k } at the output of the AID the LMS algorithm processor. 

80 in the receiver 14. The SCISF 90 filters and transforms 35 The output sequence {x k } of the second pass through the 

the received sequence {y k } into an output sequence {x'*}. SCISF 90 is compared (in signal comparator 205) to the 

The received sequence {y t } is also input to a delay 305 reference sequence. As noted above, the received sequence 

and then to a secondary SCISF 300 that has the same {y k } passes through a delay 305 before being input to the 

coefficients as the primary SCISF 90 following the initial SCISF 90 for the second pass. The delay 305 compensates 

training. Such a configuration allows periodic training to be 40 for the processing delay in the demodulation/decoding chain 

performed without disruption to the operation of the com- and the encoding/modulation chain. The comparison results 

munication system. The secondary SCISF 300 is periodi- in an enor signal e^. that is input to the LMS algorithm 

cally or continuously trained using an algorithm similar to processor 215. The training process determines coefficients 

that used for the initial training. The new coefficients of the for the SCISF 90 so that the output {x k } matches the 

secondary SCISF 300 are periodically copied to the primary 45 reference sequence {x k } as closely as possible in a least 

SCISF 90. squares sense, i.e., the mean square error between the output 

In an initial training process, the output sequence {x k } of and the reference sequence is minimized. The coefficients of 

the secondary SCISF 300 would be compared (in a signal the SCISF 90 then are updated to the coefficients determined 

comparator 205) to a predetermined sequence x k stored in in the training process. 

memory. However, as discussed above, a sequence of data 50 In a second embodiment, the SCISF 90 coefficients are 

communication bits is used as a reference sequence for chosen so that the frequency response of the SCISF matches 

training rather than a predetermined sequence. As such, the a desired spectral response GJw) that seeks to minimize the 

receiver must have a way of generating a reference signal to effects of noise -bleeding and maximize system bit through - 

compare to the output of the SCISF. put. The desired spectral response G/co) is determined 

To compute the reference sequence, the receiver essen- 55 based on the signal-to-noise ratios observed in the various 

tially duplicates the encoding and modulation processes of frequency bins of the DFT 120 in the receiver, 

the transmitter using the output of the decoder 140. Because For example, an OFDM system may have M tones, N of 

the initial training has already been performed, the SCISF 90 which (m, through nv) are used. The system operates over 

output {x'jj matches the predetermined sequence {x k } a channel with analog frequency response H c (f). Referring 

closely and ISI and additive noise are minimized. Hence, the 60 again to FIG. 1, the analog noise power spectral density at 

data output of the decoder 140 closely matches the trans- the input of the receiver A/D 80 is S^f). Prior to receiver 

mitted sequence of communication data bits x*. The data bits A/D 80, the received analog signal may be filtered by an 

are input to an encoder 320, an IDFT 330, a parallel to serial anti-aliasing filter (i.e., receive filter 75) having a transfer 

converter 340, and a prefix adder 350 similar to those in the function H a (f). The effective discrete -time impulse response 

transmitter. The output sequence {x k } of this chain is input 65 (EDIR) of the transmission channel of the OFDM system 

to the LMS algorithm processor 215 and used as a reference (including the transmit filter 65 and receive filter 75) is h(n). 

sequence in the training algorithm. The output of the A/D 80 is input to a SCISF 90 having an 
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where is a constant and * denotes a convolution of the 
discrete Fourier transforms. If the noise in the bands occu- 
pied by unused tones is sufficiently attenuated by the anti- 
alias filter (receive filter 75), r|(k) is approximately: 



where x(n) is defined as: 



pi 



mj . . . m^ are the used tones, and C 3 is a constant. A vector 
of frequency magnitudes g is defined as: 



(8) 



and s^l), the 1th component of s^ is defined by: 



0) 



(10) 



10 



impulse response g(n). G(co) is the spectral response corre- 
sponding to A(n). 

The expected signal energy /*(k) observed in frequency 
bin k at the output of the DFT120, which has a length of NM, 

is: 



m N_ 



(it) 



(4) 



Since G k =\G(nm k /M)\ 2 J the minimization of the cost func- 
tion is performed over all positive values of as: 



10 



= arg min F(g) 



where C a is a constant, 1/T is the sampling frequency and 
D k is the transmitted power in frequency bin k. The 
noise power r\(k) in bin k is: 



where: 



(5) 



c?={g^.1!glhi,G^o,iii^} 



(12) 



(13) 



A variety of constrained optimization strategies may be used 
to so lve the _a boye equations for g opr 
/"""T^nee the optimll^mpulse^response-g^-and-desired^ 
20^^ctt"al^resppnse^ G^co) (which may be expressed as 
GJnmfJM) for a system having M tones) have been " 
^ ^dt^fmiHTd; ^ r training process is used to adapt the S GISF 90 - 
thateits. impulse response g matches the desired spectral 
^response. As shown in FIG. 4, the training process may be 
25 generalized as a feedback system. A reference sequence x* 

(6) is input to the system. This corresponds to inputting a 
predetermined reference bit sequence to a transmitter. The 
reference sequence passes through a transmission channel 
410 having frequency response H(f) (including the physical 

30 transmission channel and the transmit and receive filters). 
Additive noise r\ k from the transmission channel is repre- 
sented in this general model as an external input 420 to the 
system. The resulting signal y k is input to a filter 430 having 
a frequency response G(f), e.g., a SQSF. The output of the 

(7) 35 filter 430 is then passed to an adaptation processor 440, 

which computes an error signal based on the feedback loop 
450 and adapts the filter accordingly. The adaptation pro- 
cessor may, for example, use the LMS algorithm described 
above. 

40 The reference sequence x k is also input to the feedback 
loop 450, which passes the reference sequence x k through a 
scaling filter 460 with frequency characteristic Q(f). The 
frequency characteristic Q(f) of the scaling filter 460 (which 
may be expressed as a set of frequency domain scaling 
45 factors Q^) is determined so that the SQSF adapts to the 
desired spectral response. The output of the scaling filter 460 
is used a reference for the calculation of the error signal in 
the adaptation processor 440, as described above. 



The SNR in frequency bin k i k GjJs k T g for scalars r* rk and 
vectors s*. The scalars r* are defined by: 



(9) 



^ "Using the general feedback system shown in FIG." 47a 
50^SClSE~hav.ing an impulse=response-g may be trained to 
^minimize the error ||q*x-x*g*b-Ti*g|| 2 The resulting filter 
matches P 2 (co) in the frequency domain in a least -squares 
[sense, .where: 



55 



&M|//(a0| 2 +S,(w/ 



(14) 



rS^((o)" is~the~power-spectral~density~at the input of the 
eoLsystem, .H(co) is the frequency^response of the effective 
^discrete-time impulse response (EDIR) of the transmission 
' channel, S^co) is the power spectral density of the additive 
^noise, and_Q(to) is the spectral response of the scaling filter--. 
460 having impulse response q. The solution for g in the 



To determine an expression for g that maximizes system 
bit throughput, the capacity of each frequency bin k is 

approximated by log(l+SNR Jfc ). Accordingly, the optimal 65 ^equations above specifies only the magnitude of the spectral- 
spectral profile is determined by minimizing the cost func- (response-.of theJSCISF. If the SCISRis a FIR filter, a linear 
tion F, where: 4 phase- characteristic may be used. If the length of the SCISF 
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is n^, the desired values of G(w) for the frequency bins of an algorithm similar to that used for the initial training. The 

interest are: new coefficients of the secondary SCISF 300 are periodi- 
cally copied to the primary SCISF 90. 

u i haw I — m J -frMtfat ~ ] ) ) < 15 > To compute the reference sequence, the data output of the 

GdCjrMWMjAV^Wc^ )■ s decoder 140 is input to an encoder 320. The resulting 

frequency domain signal is input to a scaling filter 520 which 

The values O are defined b * applies the set of frequency domain scaling factors Q k that 

k y * causes the SCISF to adapt to the desired spectral response, 

as discussed above. The scaled signal is input to an inverse 

GA™ k /M)(S,{jnkfM)\H(jnk/M)? + (16) 10 discrete Fourier transform 330, a parallel to serial converter 

s n ijjik/M)) 340 and a cyclic prefix adder 350, resulting in a scaled 

~ s n (jnklM)H*{jnklM) reference sequence (x*q)*. The comparison of the output 

sequence {k k } to the scaled reference sequence (x*q) k 

, r _ , ^ , , . results in an error signal e* that is input to the LMS 

The values of Q t may be computed during an initial _ , ... -f- , . • * At • rt 

t . . • j j i_ • j • 11 * . j j # 15 algorithm processor 215. The training process determines 

traonmg pcnod and may be periodically updated during co K efficients F for the SIMmasty SC ISF 300 so that the output 

operation of the communication system. fAl 1( , , * / * \ 11 

As shown in FIGS. 5, 6Aand 6B, the general feedback matches the scded reference sequence (x*q), as closely 

training process may be used to perform an initial training of as P 0SSlble m a least S( l uares sense > 1 ' e - the mean *l uare 

a SCISF followed by periodic training analogous to the error between the output and the reference sequence is 

process described above with respect to FIGS. 1-3. One 20 minimized. Periodically, the coefficients of the secondary 

difference between the techniques is that a scaled reference SCISF 300 are copied to the primary SCISF 90. 

signal (x*q) A is used rather than an unsealed reference Alternatively, as shown in FIG. 6B, the periodic training 

Referring to FIG. 5, to perform the initial training, a may be performed with a single SCISF 90. In this 

predetermined sequence of bits x* is input to the transmitter. configuration, a received sequence {y k } is output by the A/D 

The transmitted signal is filtered and noise-corrupted by the 25 80 in the receiver. The SCISF 90 filters and transforms the 

transmission channel, resulting in a received sequence {y k } received sequence {y k } into an output sequence {x k }. The 

at the output of the A/D 80 in the receiver. The SCISF 90 received sequence {y^} is also input to a delay. After the 

filters and transforms the received sequence {y*} into an received sequence {y J passes through the SCISF 90, a data 

output sequence {x k } The output sequence {x k } is compared switch 360 is changed from position A to position B, 

(in a signal comparator 205) to a scaled reference sequence 30 allowing the delayed received sequence to make a second 

( x *q)*- pass through the SCISF 90. An output switch 370 also may 

The scaled reference sequence is computed from a copy be d ^ that data is nQt ou[ t durin tfae ^ 

of the predetermined sequence x k that is stored in memory s ^ addid ^ sasF mfEdsats are controUed by 

210 in the receiver. As a first step the predetermined ^ ^ t 

sequence is input to a serial to parallel converter 510 and a A c • « j ■ ^ 

DFT 515. Tne resulting^equency domain signal is -input-to 35 ^IT^T™ IT?? ^ 1 ^ nfi ? uratlon 

a scaHng talter_52fl^hich applies the set of frequency of FIG. 6 A. The data output of the decoder 140 is input to 

domain scaling factors that causes the SCISF to adapt to ™ encoder 320 « n& resulting frequency domain signal is 

Uhe<lesiredsj^c^ input to a scaling filter 520 which applies the set of fre- 

signaLis^ quency domain scaling factors Q k that causes the SCISF to 

r^llerto-serial converter340 and"~a^y^icIpreflxTadder' 40 ada P l to the desired spectral response, as discussed above. 

t 350," resulting in a scaled reference^s~eqWnce~(x*aJj^The The scaled signal is input to an inverse discrete Fourier 

^cdmparison of -the output sequence ~{x k } - to -the -scaled transform 330, a parallel to serial converter 340 and a cyclic 

reference sequence (x*q) A results in,ah erior si^al e^thatis prefix adder 350, resulting in a scaled reference sequence 

^mpuOp^th'e LMS llgorithm _processor_215-along_.with (x*q)*- The scaled reference sequence is input to the LMS 

sequence" {x k } Alternatively, a frequency dpjrMnjxference 45 algorithm processor. 

(e.g., a predetermined bit sequence that has been processed The output sequence {x k } of the second pass through the 
by a serial to parallel converter and DFT) may be stored in SCISF 90 is compared (in signal comparator 205) to the 
memory in the receiver, which would eliminate the need for reference sequence (x*q) k . As noted above, the received 
the serial to parallel converter and discrete Fourier transform sequence {y*} passes through a delay 305 before being input 
in the feedback loop. 50 to the SCISF 90 for the second pass. The delay 305 corn- 
Following the initial training, the SCISF is periodically pensates for the processing delay in the demodulation/ 
trained during operation of the communication system. A decoding chain and the encoding/modulation chain. The 
sequence of communication data bits x* is input to the comparison results in an error signal e* that is input to the 
transmitter. Referring to FIG. 6 A, the transmitted signal is LMS algorithm processor 215. The training process deter- 
filtered and noise-corrupted by the transmission channel, 55 mines coefficients for the SCISF 90 so that the output {± k } 
resulting in a received sequence {y k } at the output of the A/D matches the scaled reference sequence (x*q) k as closely as 
80 in the receiver. The SCISF 90 filters and transforms the possible in a least squares sense, i.e., the mean square error 
received sequence {y k } into an output sequence {x' k } . between the output and the reference sequence is minimized. 

The received sequence {y k } is also input to a delay 305 The coefficients of the SCISF 90 then are updated to the 

and then to a secondary SCISF 300 that has the same 60 coefficients jletermin ed in the training process, 

coefficients as the primary SCISF 90 following the initial [^Sj^third embodimlEt^^ 

training. The secondary SCISF 300 provides output ^'Ie1^fc6'& 

sequence {x k }, which is compared to a reference sequence \ for a communication channel having a particular noise 

during the periodic training process. Such a configuration ^profile. As discussed above, a CP is added to each symbol 

allows periodic training to be performed without disruption 65 ^prior to transmission through the communication channel to 

to the operation of the communication system. The second- ''reduce the effects of ISI. However, because the CP consti- 

ary SCISF 300 is periodically or continuously trained using Ltutes redundant data^ increasing the length of the CP reduces 
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^^the^noisc characteristics of the communication channel 

For a DMT communication system with M tones, the 
maximum sample rate W (samples/second) for a particular 
channel depends, in part, on the available bandwidth and 
hardware limitations. The sample rate includes communica- 
tion data and CP bits. For a CP length of n c , the maximum 
symbol rate (which includes communication data, but not 
the CP) is W/(2M+n c ). 

Before determining the optimal CP length, the SCISF 
should be initially trained to the channel. However, a com- 
munication system need not have a SCISF to employ the CP 
optimization algorithm. It is noted that the SCISF coeffi- 
cients determined during the training process do not depend 
on the CP length. The capacity of the sub-channel may be 
approximated as log(l+SNR i ) bits per second, so the number 
of bits per symbol is 2 l log(l+SNR ( ). For a CP length of n c , 
the maximum bit rate is expressed as a function of the cyclic 
prefix as: 



14 



WTn c 



10 



15 



20 



(17) 



25 



The optimal CP length is determined by computing the 
maximum bit rate for a set of candidate values of CP length 
and finding the length that maximizes B a (n c ). 

The signal to noise ratio SNR t - of each subchannel is 30 
determined by measuring the received signal and noise 
power and computing the ratio of the two. The noise power 
y i for the i* bin may be measured by transmitting a data 
communication sequence and computing the average of the 
squares of the errors measured at the output of the receiver 35 
DFT. The total received power (signal and noise) 6 t - for the 
i th bin may be measured by computing the average of the 
squares of the outputs of the receiver DFT. The signal to 
noise ratio is determined from the expression: S/y^l+SNR,-. 
Since the signal to noise ratio is determined in the receiver, 40 
the computed bit rate B^nJ must be transmitted back to the 
transmitter. The transmitter compares the bit rate to the 
values computed for other candidate CP lengths and selects 
the CP length n c with the highest maximum bit rate B £I (n c ). 

FIGS. 7A-7D and 8A-8D show performance simulations 45 
for test systems based on system parameters and test loops 
described in VDSL Alliance SDMT VDSL Draft Standard 
Proposal, Technical report, ANSI, 1998; and Very-high- 
speed digital subscriber lines: System requirements, T1E1.4/ 
97-131R1, Technical report, ANSI, 1997. The results are for 50 
a VDSL system working over test loops 2 and 6 of length 
4500 feet in the upstream direction. The system has a 
sampling frequency of 11.04 MHz. Noise is generated by 
near-end cross talk from an interfering ADSL and an inter- 
fering HDSL and white noise at a level of -140 dBm. The 55 
SCISF used in the simulations is length-15 FIR. Aversion of 
the Normalized LMS algorithm is used to train the SCISF 
during an initial training period using a predetermined 
transmitted sequence. 

FIGS. 7A-7D show the simulated system performance for 60 
a communication system having the parameters defined for 
Test Loop 2, which is 4500 feet in length. FIG. 7A shows 
channel frequency response with and without a SCISF. The 
SCISF provides a much more uniform frequency response 
across the frequency band of interest and significantly 65 
improves the signal to noise ratio (SNR) in the higher 
frequency bins. FIG. 7B is a plot of the error signal (10 



loglXfr-xJ) during the training process. The error decreases 
rapidly during the first few iterations and is nearly converged 
after only 20-30 iterations. FIG. 7C is a plot of transmitted 
power spectral density, received power spectral density and 
the additive noise power spectral density over the used 
subchannels at the output of the receiver A/D. FIG. 7D is a 
plot of SNR at the input to the receiver A/D, which is the 
maximum attainable SNR. The plot also shows the SNR at 
the output of the receiver DFT without a SCISF and the SNR 
at the outputs of the receiver DFT using an adapted SCISF. 

FIGS. 8A-8D show the simulated system performance for 
a communication system having the parameters defined for 
Test Loop 6, which is 4500 feet in length. FIG. 8A shows 
channel frequency response with and without a SCISF. FIG. 
8B is a plot of the error signal (10 loglx^-xj) during the 
training process. FIG. 8C is a plot of transmitted power 
spectral density, received power spectral density and the 
additive noise power spectral density over the used subchan- 
nels at the output of the receiver A/D. FIG. 8D is a plot of 
SNR at the input to the receiver A/D. The plot also shows the 
SNR at the output of the receiver DFT without a SCISF and 
the SNR at the outputs of the receiver DFT using an adapted 
SCISF. 

Other embodiments are within the scope of the following 
claims. 
What is claimed is: 

1. A method for adapting a primary impulse shortening 
filter in a multiple carrier communication system, the 
method comprising: 

providing a secondary impulse shortening filter, 
comparing an output signal of the secondary impulse 

shortening filter to a reference signal to compute an 

error signal, 

computing coefficients of the secondary impulse shorten- 
ing filter in an adaptive processor based on the error 
signal, and 

replacing coefficients of the primary impulse shortening 
filter with coefficients of the secondary impulse short- 
ening filter. 

2. The method of claim 1, further comprising: 
decoding an output signal of the primary impulse short- 
ening filter to form output data, and 

encoding the output data to form the reference signal. 

3. The method of claim 2, further comprising applying a 
discrete Fourier transform to the output signal of the primary 
impulse shortening filter prior to decoding the output signal. 

4. The method of claim 2, further comprising applying an 
inverse discrete Fourier transform to the encoded output data 
in forming the reference signal. 

5. The method of claim 1, further comprising: 
receiving a digital signal from an output of an analog to 

digital converter, 
applying the digital signal to the primary impulse short- 
ening filter, 

delaying the digital s,ignal to produce a delayed digital 
signal, 

applying the delayed digital signal to the secondary 
impulse shortening filter, and 

applying the delayed digital signal to the adaptive pro- 
cessor. 

6. The method of claim 2, further comprising scaling the 
encoded output data with a set of scaling factors in forming 
the reference signal. 

7. The method of claim 6, wherein the scaling factors are 
determined by: 
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measuring received noisejpower spectral density, 
computing a desired spectral response based on the mea- 
sured noise power, and 
computing the scaling factors so that the coefficients 
computed in the adaptive processor provide the sec- 5 
oodary impulse shortening filter with a spectral 
response that matches the desired spectral response. 

8. The method of claim 7, further comprising applying a 
discrete Fourier transform to the output signal of the,primary 
impulse shortening filter prior to decoding the output signal 30 
and measuring the noise power spectral density at an output 

of the discrete Fourier transform. 

9. The method of claim 6, further comprising applying an 
inverse discrete Fourier transform to the scaled, encoded 
output data. 35 

10. A method for adapting an impulse shortening filter in 
a multiple carrier communication system having a spectrally 
constrained impulse shortening filter, the method compris- 
ing: 

scaling predetermined reference values with a set of 

scaling factors to form scaled values, 
applying an inverse discrete Fourier transform to the 

scaled values to form a reference signal, 
comparing an output signal of the spectrally constrained 25 

impulse shortening filter to the reference signal to 

compute an error signal, and 
computing coefficients of the spectrally constrained 

impulse shortening filter in an adaptive processor based 

on the error signal. 30 

11. The method of claim 10, wherein the scaling factors 
are determined by: 

measuring received noise power spectral density, 
computing a desired spectral response based on the mea- 
sured noise power, and 35 
computing the scaling factors so that the coefficients 
computed in the adaptive processor provide a second- 
ary impulse shortening filter with the desired spectral 
response. 

12. The method of claim 10 further comprising: 40 
receiving a data signal from an output of an analog to 

digital converter, applying the data signal to the spec- 
trally constrained impulse shortening filter, and apply- 
ing the data signal to the adaptive processor. 

13. A method for adapting an impulse shortening filter in 45 
a multiple carrier communication system having a spectrally 
constrained impulse shortening filter, the method compris- 
ing: 

receiving a digital signal, 

applying the digital signal to the spectrally constrained 50 

impulse shortening filter to obtain a first output signal, 
computing a reference signal from the first output signal, 
delaying the digital signal to produce a delayed digital 
signal, 

applying the delayed digital signal to the spectrally con- 
strained impulse shortening filter to obtain a second 
output signal, 

comparing the second output signal to the reference signal 
to compute an error signal, and 60 

computing coefficients of the spectrally constrained 
impulse shortening filter in an adaptive processor based 
on the error signal. 

14. The method of claim 13, wherein computing the 
reference signal comprises: 6S 

decoding the first output signal to form output data, and 
encoding the output data to form the reference signal. 



55 



15. The method of claim 14, wherein computing the 
reference signal further comprises applying a discrete Fou- 
rier transform to the first output signal prior to decoding the 
first output signal. 

16. The method of claim 14, wherein computing the 
reference signal further comprises applying an inverse dis- 
crete Fourier transform to the encoded output data, 

17. The method of claim 14, wherein computing the 
reference signal further comprises scaling the encoded out- 
put data with a set of scaling factors. 

18. The method of claim 17, wherein the scaling factors 
are determined by: 

measuring received noise power spectral density, 
computing a desired spectral response based on the mea- 
sured noise power, and 
computing the scaling factors so that the coefficients 
computed in the adaptive processor provide the spec- 
trally constrained impulse shortening filter with a spec- 
tral response that matches the desired spectral response. 

19. The method of claim 18, farther comprising applying 
a discrete Fourier transform to the first output signal prior to 
decoding the first output signal and measuring the noise 
power spectral density at an output of the discrete Fourier 
transform. 

20. The method of claim 17, further comprising applying 
an inverse discrete Fourier transform to the scaled, encoded 
output data. 

21. A multiple carrier communication system comprising: 
a primary impulse shortening filter connected to receive a 

digital signal and to accept coefficients, 
a secondary impulse shortening filter connected to receive 
the digital signal, output an output signal, and pass 
coefficients to the primary impulse shortening filter, 
a reference signal generator configured to output a refer- 
ence signal, 

a comparator connected to compare the output signal and 
the reference signal and output a resulting error signal, 
and 

an adaptive processor that computes coefficients for the 
secondary impulse shortening filter based on the error 
signal. 

22. The multiple carrier communication system of claim 

21, further comprising: 

a discrete Fourier transform connected to receive an 
output signal of the primary impulse shortening filter, 
and 

a decoder connected to receive the transformed output 
signal from the discrete Fourier transform. 

23. The multiple carrier communication system of claim 

22, wherein the reference signal generator comprises an 
encoder connected to receive output data from the decoder. 

24. The multiple carrier communication system of claim 

23, wherein the reference signal generator further comprises 
a scaling filter connected to scale the output data from the 
encoder using a set of scaling factors. 

25. The multiple carrier communication system of claim 

24, wherein the scaling factors are determined by: 
measuring received noise power spectral density, 
computing a desired spectral response based on the mea- 
sured noise power, 

computing the scaling factors so that the coefficients 
computed in the adaptive processor provide the sec- 
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ondary impulse shortening filter with the desired spec- 
tral response. 

26. The multiple carrier communication system of claim 
24, further comprising an inverse discrete Fourier transform 
connected to receive the scaled output signal from the 5 
scaling filter. 

27. A multiple carrier communication system comprising: 
a spectrally constrained impulse shortening filter con- 
nected to receive a digital signal and to accept 
coefficients, 10 

a reference signal generator configured to output a refer- 
ence signal, 

a memory to store the reference signal as a predetermined 
signal, 15 

a discrete Fourier transform connected to receive the 
reference signal from the memory, 

a scaling filter connected to scale the reference signal 
using a set of scaling factors, 
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an inverse discrete Fourier transform connected to receive 
the scaled reference signal, 

a comparator connected to compare the digital signal and 
the scaled reference signal and to output a resulting 
error signal, and an adaptive processor that computes 
coefficients for the spectrally constrained impulse 
shortening filter based on the error signal. 

28. The multiple carrier communication system of claim 
27, wherein the scaling factors are determined by: 

measuring received noise power spectral density, 

computing a desired spectral response based on the mea- 
sured noise power, 

computing the scaling factors so that the coefficients 
computed in the adaptive processor provide the spec- 
trally constrained impulse shortening filter with the 
desired spectral response. 

* + * + * 
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